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COMBUSTION SCALING AND MODELING GAS MIXING DATA IN 
NRL 5000-LITER FACILITY WITH FLOW OBSTACLE 


INTRODUCTION 


Carhart and Fielding [1] proposed a fire suppression technique of inerting pressurizable spaces by 
nitrogen pressurization. This technique is attractive since an atmosphere is created that suppresses fire 
yet supports life [2]. Other advantages include the inertness of nitrogen to electrical wiring, machines, 
and air purification devices as well as, in the use of submerged submarines, the relative ease of return- 
ing the atmosphere back to normal. There are two requirements for successful utilization of this tech- 
nique. First, the fire must be suppressed quickly to avoid excessive production of toxic gases and pres- 
sure due to combustion and heat production. Second, the addition of nitrogen must occur with 
sufficient turbulence to provide adequate mixing of the pressurant and resident gases. The feasibility of 
nitrogen pressurization has been demonstrated in suppression studies of liquid pool fires in 270- and 
5000-liter chambers [3-5]. To examine possible stratification of resident and pressurant gases during 
rapid pressurization, concentration histories of pressurant gas are inferred at specified points in the 
chamber space. Early studies without combustion and flow restriction to gases have been performed at 
the NRL 5000-liter facility by Stone et al. [6-11] and in a 1/6-scale (geometric) model at the University 
of Washington by Professor R. C. Corlett’s group (NRL Contract N00014-75-C-0185). Our present 
aims are: to study rates of mixing of the pressurant gas with the resident gas during pressurization; to 
discover how various parameters such as geometry, obstacle type and position, nozzle design and loca- 
tion, and pressurization time affect these rates; and to seek scaling rules which predict concentration 
history profiles in a geometrically scaled 324 cubic meter chamber, about one-half scale to most subma- 
rine compartments. Presented in this report are gas-mixing data in their reduced form from the NRL 
5000-liter mixing chamber where an obstacle is present to impede gas flow. Fine wire thermocouples 
are used to measure gas temperatures from which pressurant concentration histories and thus mixing 
rates are inferred [12-14]. 


EXPERIMENT AND PROCEDURE 
Apparatus 


Figure 1 schematically shows the experimental arrangements. Tank 1 is a steel, cylindrical pres- 
sure tank with two convex ends, and it lies horizontally. Its volume is 2260 liters with a length of 2.33 
m, a diameter of 1.07 m, and a mean-wall thickness of 9.5 mm. Tank 2 is a steel, cylindrical modified 
decompression chamber. Its volume is 5000 liters with a length of 2.74 m, a diameter of 1.55 m, and a 
mean-wall thickness of 11.1 mm. A hatch on one end opens inwardly, allowing access to the interior. 
An air lock is centered on the opposite end. The chamber was modified to allow numerous gas- 
collecting tubes, wires, thermocouples, and pipes to penetrate the walls. The original interior nozzle 
system, lighting system, and 6.4 mm steel floor plates remain; all other hardware has been removed. In 
addition, two thin aluminum plates have been laid on the floor plates to mount the obstacle and alumi- 
num rods which position Type K, chrome-alumel thermocouples at desired locations within the 
chamber (see Fig. 2). Pressure transducers are connected to both tanks to indicate pressure during the 
experiments. A more detailed description of the equipment can be found in NRL Memorandum 
Reports by Stone et al. [6-11]. 


Manuscript submitted September 8, 1981. 
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Fig. 1 — Schematic of apparatus 
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Fig. 2 — Plan and end views of Tank 2, showing cylindrical coordinates 
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In this work, an additional, solid, rectangular obstacle is introduced to impede gas flow and simu- 
late clutter. The obstacle is stainless steel, 0.460 m square top and 1.0 m tall. Thin rectangular wings 
protruding from the sides at the bottom enable the obstacle to be secured to the floor plates. Figure 3 
shows locations of the obstacle, the nozzle, and the thermocouples. 


KEY: 
x THERMOCOUPLE POSITION 
NOZZLE OPENING 


Fig. 3 — Interior view of Tank 2 


For convenience, a cylindrical coordinate system (r,@,z) was chosen to indicate chamber interior 
positions in Tank 2. As Fig. 2 shows, a vertical plane through the center of the cylindrical tank normal 
to its center axis is taken as the z = 0 plane, with the air lock selected as the positive z direction and ; 
the hatch as the negative z direction. Obstacle, thermocouple, and nozzle arrays are characterized by aoe 


the coordinate system. a ‘4 
ww 

Table 1 summarizes the test sets in terms of experimental variables. Configurations 1 to 3 vary ae 

only in respect to their nozzle arrays, and consequently pressurization times. Filling pressures, obstacle Le 
location and type, and in situ temperature measuring points are left as a control. Table 2 details the ‘. Be 
descriptions of nozzle and thermocouple arrays and obstacle locations. a 
3: 

Procedure _ 
The procedures for all test runs were identical. Tank 1 is charged with clean, dry air to a pressure ‘. 

of 689 kPa (6.8 atm); the air temperature is allowed to equilibrate for 5 min with the tank walls. Then, ae 
with Tank | pressure at about 657 kPa (6.5 atm) (ambient temperature) and Tank 2 pressure at 101.3 oa 
kPa (1 atm) (ambient temperature) both Doric data logging systems are started by remote switch. The ne 
data loggers are allowed to run for approximately 30 s. Simultaneously, the actuated dump valve es 
between Tanks 1 and 2 is opened and an event mark is recorded on a channel on a Doric system. Tank ae 
2 is pressurized to 202.6 kPa (2 atm) absolute. Then the valve is closed and another event mark sig- ort 
nals the termination of the gas dump. The Doric system records the mixing processes for 15 to 20 Dye 
additional minutes and is then deactivated. ° 
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Table 1 — Summary of Test Sets 
[Test Configuration _{ 1 | 2 | 3 | 
Nozzle Array I III 


Nozzle Diameter (mm) 


Thermocouple Array 
Obstacle Location 


Filling Pressures 


Mean Dump Time (s) 


Number of Good 
Experiment Runs 


Run Numbers 
Data Assigned 


oa 
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a ae 
elt at te 


Table 2 — Nozzle and Thermocouple Locations : 
Nozzle Arrays: # Nozzle (r(m), @(rad), z(m)) s 
1 (1)_—(0.483, 7/2, 0) | 
(2), (0.483, + /2, 0) eo: 
(9.483, 2/2, —0.914) ri 
mm = (3)-—«(0.483, /2, 0.914) 4 
(0.483, 2/2, 0) z 
(0.483, /2, -0.914 7 
Obstacle Location: (r(m), @ (rad), z(m)) 4 
(Obstacle coordinates on horizontal diameter plane) iy 
1 (0.076, 0, 0.229), (0.533, 0, 0.229), (0.076, 0, 0.686), (0.533, 0, 0.686 : | 
Thermocouple Array: (r(m), @ (rad), z(m)) pa 
I (1) (0.076, 0, 0.991) ee 
(2) (0.152, 0, 0.991) al 

(3) (0.229, 0, 0.991) i] 
{4) (0.229, 0, 0.914) betas 
(5) (0.229, 0, 0.838) _ 
(6) (0.229, 0, 0.762) a 
(7) (0.152, 0, 0.762) or 
(8) (0.076, 0, 0.762) ae 
(9) (0.076, 0, 1.152) - 8 7 
(10) (0.152, 0, 0.152) aos 
(11) (0.229, 0, 0.152) = 
(12) (0.229, 0, 0.076) os 
(13) (0.229, 0, 0) 25 
(14) (0.229, 0, -0.076) : 
(15) (0.152, 0, —0.076) a 
(16) (0.076, 0, -0.076) » 
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The data recorded on the two Doric system tapes in EBCIDIC character set are converted to 
ASCII and merged together programmatically on a Hewlett-Packard 21MX minicomputer. A set of pro- 
grams [15] are then used to reduce and plot the data. 


DATA ANALYSIS 


Corlett et al. [14] have provided a method of transient thermodynamic analysis of the temperature 
data to infer statistical scale-modeled pressurant mixing at each thermocouple site. A brief synopsis of 
their analysis is provided here, but a complete description can be found elsewhere [12,14]. The desired 
results are a history of local pressurant mole fraction, X,, at each interior thermocouple location in the 
chamber during a scaling run. 


Two phenomena complicate this analysis. First, the chamber contents are far from adiabatic and 
heat transfer from the gas to the walls is dependent on position. To account for this heat transfer and 
still avoid having to quantify position dependence, it is assumed that the molar gas-to-wall heat transfer 
coefficient is independent of interior position. The validity of this assumption seems allowable for 
strongly turbulent flow regimes and depends on how well the regime averages along the inside surface 
of the chamber. Second, the pressurant temperature is time dependent after entering the chamber. 
While the temperature of the pressurant gas increases as the chamber is being filled, pressurant gas 
dumped towards the end of the filling time will tend towards a higher temperature than its predecessing 
pressurant gas. Also, the temperature of a parcel of pressurant increases relative to its value at injec- 
tion during compression. It is assumed only that the pressurant age distribution is the same at every 
interior point only to avoid the otherwise necessary quantification of time dependent flow regime. This 
assumption is not too far in error for times not too near initial valve opening and after mixing is com- 
plete. 


The local heat transfer rate per mole of gas, 0, is 
o- [28|(r- 1 W) 


where R is the gas constant, T is the average temperature in the enclosure, and 7p is the initial 
temperature. Since turbulent convective heat transfer is nearly proportional to flow velocity which in 
turn is roughly proportional to dinN/dt during filling, we define a characteristic time @ during pressuri- 
zation (0 < t+ < 1¢,) and after pressurization (t > ¢,), where ris anytime, ¢, is the enclosure filling time, 
and @. is the characteristic time taken just prior to the completion of dumping. We define: 


-1 
og |2au| for0<1<t, (2) 
and 2a 
6.4 [aa fort > t, (3) 
= dt t—t,>¢ 


where N is the number of moles of gas in the chamber. The dimensionless parameter 8 is used to 
characterize a molar heat transfer coefficient so that Eq. (1) is satisfied. A control volume energy bal- 
ance shows that during compression 


yT, — T — (aT/dinN) 


‘ (y - IMT - To) 


Octet (4) 
and during postpressurization 


_ din(T ~ T))/dt 2 
(y — 1)/0, oe 


where y is the specific heat ratio of air (C,/C,) and 7, is the pressurant supply temperature. 


(5) 
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The mixing of gases which results in a parcel of gas at temperature 7; at any point and time can 
be determined as if a parcel of resident air at temperature 7, and a parcel of pressurant gas at tempera- 
ture 7, mix instantaneously and adiabatically at constant pressure. The fraction of pressurant gas in the 
mixed parcel is X;. The assumption of the pressurant age distribution being the same at every point 
allows an internal energy summation over the entire chamber 


NT =(N-~N))T, + Nola (6) 
which results in 
yo BoB [- BF) [3-8 a) 
T, — T, Nj\| 7,-T T,-T 


where X is the molar average pressurant fraction for the entire chamber contents. The average tem- 
perature of pressurant gas, 7,, is determined by a numerical integration over the chamber contents as 
shown 


aT, = a T,(dinP) — pl(T, — T/y)dinP — T,dinT] (8) 


where zw = T, ~— To/ T- Ty and T, is the supply temperature of pressurant measured at a nozzle exit. 
Program CINFR[15] carries out this integration using the following initial boundary conditions: 
T, = T = Tp), T, = Ty, P = Po, and w= wo = 1. Thus the pressurant mole fraction at any pvint 
inside the chamber can be calculated for times both during and after pressurizaiton from only a tem- 
perature profile and measured pressure data. 


A scale-modeling hypothesis is also produced. It states that, for geometrically similar enclosures, 
a dimensionless pressurant deviation, é, is a unique function of dimensionless time, 7, at homologous 
model and prototype points regardless of enclosure pressure level or pressurization rate [14]. Dimen- 
sionless time is defined for the filling process as 


r= + O<t<% (9) 
and for postfilling as F 5 ‘é 
T= 0. +1 ¢>4 (10) 
where @, = —[In(1 — X)/ dt];*, can be used. The dimensionless pressurant deviation is defined as 
X-X 
[= ; (11) 


c 

where X is the inferred pressurant mole fraction at any time 7; X is the molar average pressurant frac- 
tion where all pressurant added to the chamber in time + is perfectly mixed with the chamber resident 
gas; and X, is the maximum, final value of X at time ¢ = 4. Notice the pressurant deviation is positive 
when the chamber gas is pressurant rich, negative when pressurant poor, and zero when completely 
mixed. 


This report provides a relationship for pressurant gas concentrations as a function of time and 
position so that: 


(1) the gaseous mixing which results from the turbulent flow around an obstacle can be 
inspected and described in some way; and 


(2) mixing can be characterized and predicted for any geometrically similar enclosure with a 
geometrically similar flow regime. 


This relationship is statistically averaged over several scaling runs for each test flow configuration. 
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RESULTS 


As mentioned earlier, a minicomputer is utilized to systematically reduce and analyze the data for 
each test configuration. Table 3 indexes all data charts produced by computer software. A first pair of 
programs is responsible for reducing the data in each run into tabular form. (See Tables 4, 10, and 16 
for a summary of raw data.) The first column is time in seconds where ¢t = 0 is the time when pressuri- 
zation of Tank 2 begins. Notice that the pressure and temperature of the pressurant gas in Tank 1 is 
only noted during pressurization. The pressure in Tank 2 (noted as Fill tank) is given during the entire 
run with the in situ thermocouple. Although the data loggers collect almost three compete scaus of 
data per second, data points are averaged down to integer-second points. This was necessary since each 
data logger scans at its own intrinsic speed although their internal clocks are initially synchronized. All 
temperatures are reported in degrees Kelvin and pressures in megabars. A table showing the location 
of each thermocouple in cylindrical coordinates is provided on the reduced data tables for convenience. 


Table 3 — Index of Data Tables 


figuration N 
dae aaah 


Table Number Table Heading 
10(A-L) 16(A-I) 


Temperature Data Tables 
(per scaling run) 


11(A-D 17(A-H) | Inferred Pressurant 
Distribution (per 


scaling run) 


Mean Value of All 
Quantities (per scaling 
configuration) 


Normalized Mean Local 
Pressurant Fractions 


Standard Deviation of Mean 
Local Pressurant Fraction 


Deviations of Mean Local 
Pressurant Fractions 


Note: Tables 4-21 appear at end of report. 


A second pair of programs run on the minicomputer manipulates the reduced data tables from 
one test configuration into two new sets of tables. In the first set, for each test run, a table of inferred 
pressurant distribution is calculated. Time is given in seconds, temperature is in degrees centigrade, 
and other quantities are dimensionless. The first six columns give: time 4 molar mean temperature 7, 
mean resident air temperature 7,; mean pressurant air temperature 7,;, the dimensionless coefficient 8; 
and the ratio 8/@. Columns 8 through 20 are the calculated pressurant_fractions X; at each of the 13 
respective thermocouples, and column 7 is the mean pressurant fraction X (see Tables 5, 11, and 17). 


The second set consists of four groups of statistical tables calculated for each single test 
configuration. The first group of statistical tables (6, 12, and 18) consists of the mean values of all 
quantities described in the previous tables. A last column at dimensionless time, 7, is listed in column 
21. The second group of statistical tables (7, 13, and 19) lists the normalized mean local pressurant 
fractions for each of the 13 thermocouple locations normalized by the mean pressurant fraction at the 
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end of dumping (1 = ¢,). Other tables (8, 14, and 20) indicate the calculated values for the statistical 
standard deviation of the mean local pressurant fractions. This value might be used to indicate a 
magnitude of precision in the preceding statistical charts. A last group of tables (9, 15, and 21) displays 
the pressurant deviation é, for each normalized time 7, in dimensionless increments of +0.05. 


DISCUSSION 


The reduced data for three obstructed air mixing configurations have been presented and can now 
be discussed to some degree. While it is not possible to verify the two modeling assumptions directly 
from the results, inspection of the Inferred Pressurant Distribution charts shows an almost constant 
value of beta during the filling interval for each run. This observation 1: in accord with the expected 
heat transfer between the walls and the surroundings for highly turbulent convective flow. While a 
convective heat transfer coefficient is not known for any position along the chamber wall, an averaged 
value over the entire chamber surface is presented and is now known to be constant while its value is 
required [15]. 


The thermodynamic treatment of the data no longer describes reality after complete mixing has 
been achieved. Now, variations in the temperature profile may be caused by several influences. First, 
flow inside the chamber eventually slows into a random laminar character as the gas kinetic energy is 
eventually dissipated by viscous drag forces. Second, the chamber itself may act as a heat sink for its 
contents and cause the gaseous contents to begin to cool. Last and least of consequence, the datalog- 
ging system is known to drift (<0.6°C) in an isothermal environment. Therefore, at some time after 
mixing is complete, a mixing profile is expected to show a deterioration of perfect mixing. 


To characterize a mixing profile, we have plotted a dimensionless mixing deviation, 5, as a func- 
tion of time. 


6, = = == 7. (12) 


This value is effectively the difference between the actual pressurant concentration and its theoretical 
value if perfectly mixed normalized by its perfectly mixed value. Notice it is positive when pressure 
rich, negative when pressurant poor, and zero when theoretically mixed. Figures 4 to 6 show this plot 
at four locations for the test configurations 1 to 3, respectively. These graphs are the result of statisti- 
cally averaged values of mixing deviation for runs with identical dump times for a given configuration. 
Inspection of their corresponding statistical standard deviation shows that the assorted peaks and curves 
before complete mixing are meaningful with good precision. For the purposes of the plotting software 
[16], it was necessary to indicate the right-hand limit of the commence of value opening as zero time in 
deference to previous practice. 


Test configuration one is clearly shown to have the slowest, most stratified mixing profile. Here, 
one centered nozzle probably forms a wake as the discharged gases are reflected from the chamber 
floor. Flow vizualization films by Stone show that the pressurant gas ceases to jet almost a foot above 
the floor underneath the nozzle. During the first few seconds, a thermocouple’s mixing curve initially 
agrees with intuition, i.e., thermocouples close to the nozzle are pressurant rich and those further away 
are pressurant poor. As a wake develops and scatters, the chamber contents become temporarily 
stratified. The pressurant concentration at a location may alternate being pressurant rich and poor. 
Thermocouple number five is shown to be located inside the wake’s origin. Since it is close to the noz- 
zle it remains pressurant rich much longer than the other locations. 


Configurations two and three exhibit a different approach to complete mixing. In these cases, the 
flow circulation from more than one nozzle appears to become so turbulent that a pattern much more 
complicated than a wake may exist and mixing is much more rapid with respect to their dumping times. 
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Fig. 4 — Test configuration one: one centered nozzle. 
Dump time = 33.0 
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Fig. 5 — Test configuration two: two centered nozzles. 
Dump time = 17.0 
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Fig. 6 — Test configuration three: three centered nozzles. 
Dump time = 12.0 


A scale-modeling hypothesis states for geometrically similar enclosures, the dimensionless pres- te 
surant deviation is a unique function of dimensionless time at homologous model and prototype points mG 
regardless of pressure level or pressurization rate [14]. These quantities have been defined earlier and 
are plotted in Figs. 7 to 9 corresponding to test configurations 1, 2, and 3. By hypothesis, these curves 
should look the same for any size chamber with similar geometry. In these graphs, the left-hand limit a 
of the commence of valve opening is taken as zero dimensionless time, tau. Tau is unity when the ges 
valve has been fully closed, i.e., one dump time. The data for these plots come from the same scaling z 
runs as the previous graphs after being statistically combined by program CINST. Mixing profiles of 
the same type are seen in these graphs; however, the mixing which was done during valve opening is aa 
now observable. For our most hindered location studied, our data show that mixing is complete within eure 
7 = 1.1. A more detailed description of mixing field would require more than the thirteen thermocou- ; i‘ 
ple locations studied. Es e 


Since these configurations vary by the effects of nozzle number and location, the flow and mixing 
profiles vary. The degree of clutter by which a particular interior location is characterized is seen to 
have two effects. One effect is how well enclosed, or isolated a location is from the pressurant flow. 
This effect hinders mixing. A second effect is observed where gas is forced around proximal obstacles. 
Here, clutter may aid mixing as eddies are formed inside flow boundary layers or outside jet gas 
entrainment boundary layers. This might be considered similar to baffling in a mixing tank. A method 
of quantifying degree of clutter which incorporates both of these effects is presently intuitive, at best. 
This method of mixing analysis is instrumental in measuring degree of clutter and final mixing time. 
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CONCLUDING REMARKS 


Presented are raw and reduced data from three cluttered air mixing experiments. The measure- 
ment of pressurant concentration distribution is based on temperature measurements in the mixing 
chamber. Although this is a relatively inexpensive method for gathering data from a number of enclo- 
sure points, it suffers from (a) the loss of precision due to modeling and thermodynamic assumptions 
and (b) the requirement for extended computer time to reduce the raw data. The first drawback can be 
alleviated by employing a more accurate method of monitoring pressurant inflow rate, and the second 
has been worked out as computer programs have become more versatile and efficient. 
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